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oi THEORETICAL ANALYSIS OF THE LATERAL STABILITY 

L ? OF A GLIDER TOWED BY TWIN PARALLEL TOWLINES 

By Marvin Pitkin and Marion 0. McKinney, Jr. 

SUMMARY 



A theoretical analysis of the lateral stability char- 
acteristics of a gltdftC towed "by two parallel towline* 
has "been made and correlated with previously reported dy- 
namic flight tests of a glider model in the NACA free- 
flight tunnel, Caiculatiors were made for a range of ef- 
fective dihedral angles frcm -6° to 14°, a ranee of tow- 
line lengths from one to three glider spans, lift cceffi- 
cient? of 0 6 30 and 0,75, and a range of glider lift-drag 
ratios from 4,5 to 21,0, The results of the theoretical 
analysis indicated that a glider towed "by twin parallel 
cables possesses sufficient stability for satisfactory 
pilotless towed flight. 

The stability of the glider was found to he primarily 
dependent upon the dihedral angle of the glider and the 
length of the towlines. Both unstable oscillations and 
divergences were found to occur either at negative or at 
large positive effective dihedral angles. A maximum of 
stability was indicated at moderate positive dihedral an- 
gles, -he stability of the glider war decreased as the 
towline length and lift coefficient increased. Large 
sidewise displacements of the glider relative to the tug 
reduced the stability, particularly when the dihedral an- 
gle was large. Agreement between theoretical and flight 
results was very satisfactory. 

INTRODUCTION 



Load-carrying gliders may be used to carry troops or 
cargo to supplement the carrying capacity of existing air- 
planes without seriously decreasing the efficiency of the 
tug. In order to realize full benefit from such an ar- 
rangement it is desirable to have inherent stability in 
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the glider system. An adequate decree of stability would 
eliminate the need for a pilot in the glider or would at 
least relieve the pilot of the necessity of giving con- 
stant attention to the control?. Although some success 
has "been obtained with a single towline glider system, 
considerable difficulty has been experienced in obtaining 
a sufficient degree of lateral stability. A somewhat 
simoler system for eeourihg stability appears to be one 
involving twin parallel towlines attached near the wing 
tips of the glider to restrain the glider from yawing. 
A preliminary investigation (reference l) in the VAC A ■ 
free-flight tunnel indicated that a stable pilotless towed 
glider system was possible with twin parallel towlines. 

In order to obtain an understanding of the operation 
of this, yawing restraint in providing lateral stability 
and to allow an extension of the results to other glider 
configurations, a theoretical analysis has been made of 
the lateral stability of the system. Calculations "ere 
made for the Bristol % pw~€arget glider ,, Sk3et ,! to corre- 
late the theory with the tests of reference 1. 



SYMBOLS 



* weight of the glider, pounds 

IB mass of the glider, slugs 

b span of the glider, feet 

*b i distance along Y- axis from the center of gravity 

of the glider to the towline attachment 
points, feet 

3 wing area of the glider, souare feet 

X,Y,Z cent roldal axes of glider 

a length of the towlines, feet 

a 1 length of towlines in glider span lengths (a/h) 

radius of gyration about the X axis, feet 



A ' ,B\C« 



space axes parallel to the body axes of the tug 
intersecting midway between the towline at- 
tachment points on the tug 
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£ angle "between t owline and A 1 -axis, radians 

V true airspeed , fee': per second 

W 

T~) 3 angle of sideslip 9 radians 

t 

y sidevi s e displacement of glider center of gravity 
along the 3 '-axis, feet 

y 1 sidewise displacement in glider span lengths (y /t) ) 

v sideslip velocity, feet per second 

<t> angle of tank, radians 

p rolling angular velocity f radians ver second 

\|/ angle of yaw, degrees 

D drag of glider , pounds j differential operator ( — ) 

\ds/ 

T t owl in? tension, pounds 



C L lift coefficient 



/ 



ift\ 
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/drag\ 

Cp) drag coefficient / — j 

\|sv J 

p density of air, slugs per cu"bic foot 
t time , seconds 

n airplane relative— density factor (m/pSb) 
L rolling moment of glider, foot-pounds 
Y lateral force of glider, pounds 

Cj rolling-moment coefficient i - — \ 

\3 
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lateral-force coefficient / - \ 

\2 / 
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v 

x 

V 



rate of change of lateral force with sidewise 
vel ocity (dY/dv ) 



Y\j/ rate of change of lateral force With angle of 

yaw (dY/c40 

(Yy) 1 rate of change of lateral force irduced by tow- 
line tension with glider sidewise displace- 
ment (dY/dy^ 

(^fy)? rate of change of lateral force- due to yaw pro- 
duced "by sideslip and roll with angle of 
roll CdT/M) 9 

Lp rate of change of rolling moment with rolling 

v e I o c i t y ( o L / 6 p ) 

rate of change of rolling moment with sideslip 
velocity (<3L/6v) 

Z»^h rate of change of rolling moment with angle of 

yaw (dL/dty) 

(I»<J>) i rate of change of rolling moment induced "by tow- 
line tension with angle of roll (dL/dO)} 

(l»(t))a rate of change of rolling moment due to yaw pro- 
duced by sideslir> and roll with angle of 
roll (dL/d<t>>3 

(lk{))a effective rate of change of rolling moment due to 
yaw produced by sideslip and roll with angle 

of roll (*£) 

(Y(h) 3 effective rate of change of lateral force due to 
yaw produced by sideslip and roll with angle 

of roll f ~ , 
\6<t> 1 

(Cyg) rate of change of lateral-force coefficient with 
angle of sideslip, per radian (dCy/dP) 

(Cy y I )i rate: of change of lateral-force coefficient in- 
duced by towline tension with sidewise dis- 
placement (dCy/dy ? 
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( C y<£ ) effective rate of change of lateral-force coef- 

ficient due to angle of yaw produced "by 
sideslip and roll with angle of roll 

~1 rate of change of rolling-moment coefficient 

with angle of sideslip, per radian (dCi/d&) 

coef f icient 



coef f i ci ent 

th anfrle of 



^l&'a effective rate of change of rolling-moment co- 

efficient due to angle of yaw produced "by 

sideslip and roll with angle of roll ^InT J 

s forward distance traveled in s"Can lengths f^*t J 

\h y 

h root of the stability quart ic 

A,B,C,S,F coefficients of the stability quartic 

R South' s discriminant (BCE - S 2 - B 3 F) 

id imaginary portion of a complex root 

c real portion of a complex root or a real root 

P period of Oscillation, seconds 

t 3 / 2 time reouired for motion to damp to one-half 

arnpl i tude, seconds 

r geometric dihedral angle, degrees 

* e ff effective dihedral angle, degrees (-Cj /0.012) 



C j, rate of change of r olling-momer. t 

with helix angle \^ C y J 

(Ol^Ja rate of change of rolling-moment 

induce^ by towline tension wi 

roll ( 8 C 1 1 dO ) a 
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METHODS 



The analysis of the lateral stability of a glider 
attached to a towing airplane by twin parallel towlines 
was made "by setting up the stability equations, obtaining 
the stability derivatives when necessary, and then calcu- 
lating the stability boundaries, the period of the oscil- 
latory (periodic) mode, ard the damping of the oscillatory 
and the aperiodic components of the glider motion. 



Asetunpt 1 ons 

In order to facilitate the handling of the theoreti- 
cal investigation the following assumptions were made: 

(l) Level tow: The tug and glider were assumed to 
be in level flight with the towline? horizontal, 

(3) Straight towlines: The towlines werp assumed 
weightless and straight. This assumption is justified by 
the flight data of the investigation reported in reference 
1. Thus, no derivatives considering the weight and curva- 
ture were included. 

(3) Two degrees of lateral freedom: The glider was 
considered to be restricted by the towlines to two de- 
grees of lateral freedom. 

(4) Towpoint attachments on Y-axis: The glider was 
assumed to be attached at two points on th^ Y-axis. 



Equations of Motion 

Inasmuch as the dyadic towline system restrains the 
glider from yawing, only two degrees of lateral freedom 
need be considered in setting up the equations of r.otion 
that apply to this problem 0 The lateral-force and rolling- 
moment equatiens after a disturbance may conseauently be 
written as 

dv 5Y 
&^t* = Nfc + vt-- + towline derivatives (l) 

. »&p dL hL f N 

m&yr — = p- — + v — + towline derivatives (2) 
dt 0 "d dv 
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The problem therefore becomes or. e of evaluating the 
towline derivatives. When this step is accomplished, the 
stability equations may be set up and treated mathemati- 
cally to obtain definition of glider stability, 

T owl in e Derivatives 

Derivatives due to cable tension.- Two towline deriv- 
atives due to cable tension may be obtained from inspec- 
tion of figure 1, A sidewise displacement of the glider 
along the direction of its Y-axis is opposed by a compos 
nent of cable tension. The derivative expressing the re- 
sultant relationship is then 

(X ) mlSL) 

-Tens i on (or dr ag) 
Towline length 

m ± ■ (3) 

a 



to 
I 

1-3 



Similarly, the towline tension resists the effort of the 
glider to assume a banked attitude. The derivative de- 
fining this r elat i onship is 

Tension (or drag) x b j 
Towline length 



-Dd 3 3 



The nondimensi onal form of these equations is 



(4) 



and 



_5 f°i 



a 



(5) 



(6) 



B 



The towline derivatives 




and 




are 



the only first-order terms contributed "by the tovlines. 
Inasmuch as conventional stability theory assumes that 
only first-order terms are considered, these terms were 
the only towline derivatives to "be introduced into the 
equations of motion (l) and (2) in the early phases of 
the investigation. The glider stability was then calcu- 
lated. 

The results of these initial calculations were some- 
what in disagreement with the model flight data of refer- 
ence 1 because the calculations did not predict instabil- 
ity at the large dihed^ .1 angles at which unstable oscil- 
lations had been encountered in the flight tests. A clue 
to the correct solution was obtained from a further anal- 
ysis of data previously obtained from flights of the model 
with large dihedral angles. The model was observed to fly 
satisfactorily during the tests of reference 1 provided 
that the tunnel air-flow disturbances were small. For 
large disturbances involving large glider displacements, 
however, the model quickly developed an unstable oscilla- 
tion. This phenomenon indicated that large displacements 
were introducing unstable effects sufficiently large to 
overcome the inherent stability of the glider system for 
small oscillations. It thus appeared likely that some 
second-order effect existed which was sufficiently large 
at large values of glider displacement to induce instabil- 
ity. 

An analysis was therefore made in an effort to deter- 
mine which, if any, of the second-order terms of the 
glider system could have an appreciable effect upon glider 
stability. This analysis indicated that the second-order 
effects of the derivatives hitherto discussed would be of 
small magnitude, even for large values of glider displace- 
ment. It was observed, however, that new towline deriva- 
tives, second-order in form, arose from glider displace- 
ment and were of sufficient magnitude for large displace- 
ments to induce instability. These derivatives resulted 
from the geometric characteristics of the parallelogram 
formed by the towlines and tJie tug and glider wing spans. 
The derivation of these derivatives and the method of 
introducing them into the glider equations is discussed 
in the following section. 

Second- order derivatives due to glider position.- The 



geometric characteristics of the parallelogram formed by 
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to 



the towlines and the tug a 
if the glider sideslips an 
the towlines force the gti 
horizontal tow this action 
panying positive sidewise 
yaw introduces additional 
dependent upon the amount 
bank. If the effects of y 
lieved to "be small, are ne 
and the lateral farce t 
"by glider position are fou 



nd glider w 
d r a J> 1 ^ pel 
d e r into 

r e suit s in 
di splacemen 
stability d 
of sidewise 
awing veioc 
glectfid, th 
due to the 
nd to "be 



ings are such that, 
ative to the tug, 
yawed attitude. In 

positive yaw accom- 
t . Thi s angl e of 
erivatives that are 

displacement and 
ity, which are be- 
e-rolling moment L 
angle of yaw in-diu>e<ir~ 



L 

v 



But 



therefore , 



vii 



L = f(ct> f y)L 



(?) 



Li kevi a e , 



Y = f(0,y)Y^- 



(8) 



The exact solution for if in terms of 0 and y 
leads to a highly complex equation. An approximate for- 
mula, however, has been obtained by spherical trigonome- 
try and checked by means of a mechanical device- reproduc- 
ing the glider tew system. The eouation for angle of yaw 
is 



0 C ' tan c 



1.9 ~ C.07(a ! - 1) 



(9) 



where 



tan c = 



It should be observed that the sidewise displacement 
can never exceed the length of the towlines a f . 
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The agreement "between values of calculated "by 

means of this equation and the measured values is shown 
in figure 2, which presents sample values for 1-span tow- 
lines. The agreement is eaually good for 2- and 3-?pan 
t owl in 6 8 • 

Substitution of equation (9) in equations (7) and (8) 
gives the relationships, 



r 

1.8 - 0.07 (a f - 1 ) 



(10) 



and 



Y ss (n) 



J (a 1 ) 3 - ty 1 ) 2 \1.3 ~ 0.07 (a' * l) 

Steps must now be taken to obtain a relationship be- 
tween angle of bank 0 and the sidewise displacement of 
the glider y ! because, first, an infinite number of so- 
lutions for glider stability would otherwise exist and 
second, insertion of derivatives as a function of <t> and 
y 1 in the equations of motion would make these differen- 
tial equations nonlinear and hence unsolvable by conven- 
tional mathematical methods. 

Although the exact relation of the angle of bank to 
the sidewise displacement can be found only by a complete 
solution of the equations of motion, a satisfactory first 
approximation may be made if the towline and inertia 
forces are ignored. .The relationship of 0 to y ! may 
then be expressed by the eauation 

r>b 

ec ^ + m c h = 0 



or 



*0 ~ *2? " 4 p 



PC la = -II- Cj (12) 



The value of C? as used herein has been taken from 

p 

reference 2 as -0.4; then 
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. J% (-0.4, 



or 



03 

2 I n = & (0.2) 

Multiplying 'by V and integrating "both sides of this 
equation with respect to time yields 

yC lj3 = cpb(0.2) 



or 



Qj p = 0.2 



v ! 



ig - 0.2 (13) 



The validity of this equation was checked by means 
of measurements of values of y 1 /<J> obtained from motion- 
picture records of flight tests for various effective, di- 
hedral angles and lift coefficients. The test points 
shown in figure 3 are averages of at least seven readings 
for each condition and are in satisfactory agreement with 
calculated result s. 

Two tovline derivatives due to position may now he 
expressed as a function of either $ or y T . The calcu- 
lations in this reioort have utilised values derived as a 
f unc t i on of <$> . 

The rolling moment resulting from glider position 
may "be expressed as 

or, In nondimens i onal form, as 



•*0l 



= -f(O)0: 



(14) 
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li kewi s e , 

C Y = -\|/C Yp 

' = -f(<i>)C Yc (15) 
P 

In order to introduce the moments and forces arising from 
glider position into nondimensi onal eauations of motion, 
it is convenient to use the relations 



C 7 



= *( ci 4 ( 

It is then necessary to evaluate (\) and 

Inasmuch as C j is not a linear function of 0, 

will vary with 4>« In order to obtain an effective value 
of i til. i over any range of values of <> from an ini- 

tial value of <t> 0 to a final value of (J) , the nonlinear 
variation of C j with <t> can "be replaced by a linear 
variation which represents the work done in any cycle. 
This linear variation, effective (q l/\ , is derived as 

follows: N 

Work done in — cycle = \ — — / (<P - <t> 0 ) 

4 \ W - ®q / 



= mean 0» (<*> . # a ) (18) 

Jo 
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When <$ 0 equals 0, 

Work done in ~ cycle = C j$ 

The linear variation of C? with <t> yielding the same 
Work in ~ cycle is consequently defined as 

/_ 



Effective U l ^) = (C ^ 



2 



20 I 
<$> 

2 / C , dcj> 
s 



From equation (14) 



01 - 



Therefore , 



(is) 



(20) 



0 



where, from equation (9), 



1.8 - 0.07 (a 1 - 1) yi' 8 - y' 3 
and, from equation (13). 

The integral in eauatior. (20) may he solved in terms of 
0 or, by rearranging limit s , in terms of y j ; either 
process yields identical numerical results. Values of 

f 0 ) are expressed herein as a function of y 1 « 
\ Vg 
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Solution of equation (30) in terms of y 1 yields, when 
i s negat ive , 



•0 



5.33Cifi' r /""a a. 2 .s. | a' 



fcD 1 S Jy a . _ y « (y«%2 a O-2* 

\ C > 3 [l.S-0.07(a» -1)] y ,<; 



(21) 



When- C r is positive, the equation is of opposite sign. 
P 

The effective derivative fey/) may "be obtained in 

V / 2 • 

a similar manner or "by the relationship 

(t~\ (rf\ — (22) 

These towline derivatives are of second order and 
are dependent upon glider position. The glider stability 
must therefore he calculated for each sidewise glider po- 
sition assumed. 

The towline derivative has ^ een plotted 

against sidewise displacement for various value? of Cjg 

and is presented in figure 4 for 1-, 2-, and 3-span tow- 
lines. 



Stability Squat i on s 

Substituting the towline derivatives in equations (l) 

and (2) and rearranging terms yields the following equa- 
tion s * 

_ W0 „ vY „ y (y ) _ 0(y7) 3 = 0 (23) 

* k x 3 i? - ? l p - vL v - * ( Vi " = 0 (24) 



If the aerodynamic data in standard nondimensional 
form is to be used and if the original physical signifi- 
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I 

h3 



cance of the equations Is to "be maintained, a consistent 
nondimensi onal system is desirable. The following units 
for such a nondimensional system were considered conven- 
ient: 

Unit of length b 

Unit of mass pSh 

Unit of time b/V 

If length, mass, and time are expressed in the nondi- 
mensional form, the following nondimensional quantities 
will result : 



u = 



o 3b 
v 



s * - t 
b 



Division of the force eouation (23) by ^-SV and the mo- 
ment eouation (24) by ~SV 2 b yields the nondimensional 



equat ions 



Tcr + ( 5?7\1 - EPC Yc - P fcv t ] = 0 (25) 



k x>„3P^ p p* r/ \ / — \ 



4ul ) D ■ D- — C i - D 3C i p - 2— 

v > > ST 2V P 



= 0 (26) 



The det ermiaant of the left side of these equations yields 
the glider s t ah i 1 i ty equation of the form 

A>> 4 + BX 3 + CX 2 + EX + F = 0 (27) 

wher e 

A = 1 



—- HK*rJ ° y b - 2uCi P 



1G 



By. — y 



E = 



8(rv — y L - 

^ n / 



m 



The derivative f^-i j enters into the calculations 

\ J s 

such a way a? to cancel a term involving ( C 7 ^ . 

v <V 2 



GLIDER CHARACTERISTICS 



The dimensional and mas? charcct eristics cf the 
Bristol tow-target glider w Skeet n for which the calcula- 
tions were made are shown in the following table: 



Wing area, square feet 173.0 

Horizontal tail area, square feet 56.1 

Vertical-tail area, square feet... . . t 10.5 

Wing snan, feet 54.4 

Over-all length, feet 29.3 

Distance from center of gravity to towline- 

attachment point?, feet 15.0 

Weight, pounds 2976 

Wing loading, pounds per square foot 17.2 

Radius of gyration about X-axis, feet 5.3 

Moment of inertia about X-axis, slug-feet 2 .... 2596 



Other factors that enter into the numerical calculations 
are: 
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V 



( <k ,1 

( Oil] 



The value 



-0.40 
-0.40 



p. = 6.53 



"a- lb y 



0. 154 
Id 

Cj) = 0.0667 (at Ct = 0.3) 

Ct) ■ 0.1047 (at C L = 0.75) 

is dependent upon glider displacements 



and is presented in figure 4. 



CALCULATIONS AND RESULTS 



Inasmuch as 
t Ire s showed i C \ 



the calculation of the tovline deriva- 



and 



to "be functions of 



glider displacements, the stability was investigated over 
a range of sidewise displacements "by a step-hy-step proc- 
ess. It was first necessary to assign some limiting val- 



ue to the sidewise displac e id e n t y 



The towlines limit 



the sidewise travel to the length of the towline used "but 
inspection of the model flight records indicated that 
these maximum values were never obtained in flight even 
with the most violent gusts or divergences. A study of 
these records indicated that, in almost all cases inves- 
tigated for 1-, 2-, or 3-span towlines, the initial side- 
wise displacements of the glider center of gravity caused 
"by g^st disturbances did not exceed a value of 0 % 70 span. 
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Since the tunnel conditions simulate rather severe gust 
conditions in full-scale flight, this value was thought 
to "be an adequate limit defining the maximum initial dis- 
placement likely to occur following a gust disturbance. 

Stability Boundaries 

The theoretical requirements for dynamic stability 
of an airplane are that the coefficients of the stability 
equation and Routh r s discriminant be positive. Negative 
values of the coefficients indicate divergent aperiodic 
modes; whereas a negative value of Routh 1 s discriminant 
R (for positive coefficients) indicate* an unstable per- 
i odi c mode . 

Stability boundaries separating stable from unstable 
regions may be defined by the eauations: 

C, E, or V - 0 (aperiodic boundaries) 

R m BCE - E 3 - B 3 E - 0 (periodic boundary) 

The periodic boundary 1 = 0 was obtained by first 
evaluating the stability coefficients as defined in equa- 
tion (27) for the desired glider condition in terms of 

l 0 ] . Routh's discriminant was set up and the value of 
(^IqJ calculated for R = 0. This value was then taken 

from tne correct plot of C^'^J against y 1 (fig. 4) 

and the value of y f obtained. The values of y 1 at 
which the stability coefficients I and C became zero 
were similarly obtained. 

The aperiodic boundary S = 0 is independent of the 
pcsition derivatives and was obtained by calculating the 
value of the derivative C^g necessary to cause the coef- 
ficient E to become zero. 

The periodic boundary R ■ 0 and the aperiodic 
boundaries C = 0, E = 0, and tf = 0 are given on the 
charts of effective dihedral F eff plotted against side- 
wise displacement in figures 5, 6, and 7 for two values 
of lift coefficient and for 1-, 3-, and 3-span towlines. 
Stability ratings from flight data of reference 1 are 
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also given on these figures for different values of effec- 
tive dihedral angle. The geometric dihedral angle was 
converted to effective dihedral angle "by Leans of data 
from unpublished full-scale-tunnel tests of the glider 
made at LMAL . The full-scale data indicated that the ef- 
^ fective dihedral angle of the test glider was 1.5° lower 

than its geometric dihedral angle for a lift coefficient 
of 0.30 and 1.5° higher than its geometric dihedral angle 
for a lift coefficient of 0.75. 



Stability Roots 

The stability boundaries, although defining the var- 
ious stability regions, give no quantitative indication 
of the variation of stability within a given region. 
These characteristics may be obtained by factoring the 
stability equation for its roots X , A , X . and X, 

and using these roots to determine the period of the lat- 
eral oscillations and the time to damp to one-half ampli- 
tude for the oscillatory and aperiodic modes. 

The period in seconds of the periodic modes X 3 and 
X, was obtained by the formula 

p _ 2-n b 
d ¥ 

where id is the imaginary portion of the complex roots 

A 3 ,X 4 = 0 a, id 

The time to damp to one-half amplitude in seconds 
was calculated from the formula 

lOff*0.5 V, 

I - _ ^_ Z 

1/2 ~ C V 

0.69? b 
c v 



where c is either the real portion of the complex roots 
for the oscillatory mode or the real root X or X 

defining the aperiodic modes, 
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The reciprocal of the time to damp to one-half ampli- 
tude has "been used herein to express damping. 'This pro- 
cedure was followed "because it wa? desired to express the 
degree of stability as a direct rather than as an inverse 
function of the ordinate. Peak positive ordinates conse- 
quently indicate maximum damping and negative ordinates 
indicate negative (unstable) damping. 

The aperiodic mode \ represents a heavily damped 

convergence that is always stable and is little affected 
"by changes in glider or bowline configuration. An example 
of the damping characteristics of this mode is presented 
in figure 8 for the condition with 2-span towlines to il- 
lustrate its general n uure. The mode X is the govern- 
ing aperiodic mode and hence determines the aperiodic 
characteristics of the gl Her, 

Complete calculations of the rjeriod and damping char- 
acteristics have "been made for tow with l- f 2-, and £-span 
towlines at two values of lift coefficient and two values 
of sidewise displacement. Representative results are pre- 
sented in figures 9 to lf>« Erperimental values of the 
period of the lateral oscillations determined in the tests 
reported in reference 1 and corrected to full-scale values 
are also given on the figures for purposes of comparison. 

DISCUSSION 
Effect of Dihedral Angle 

An analysis of the stability charts (figs. 5 to 7) 
indicates that dihedral angle is an important factor in 
the determination of the stability of the glider. The 
only satisfactory stability region is that in which the 
effective dihedral angle ranges from small to moderately 
large positive value*. At angles above this region (in- 
dicated by the dashed line on charts) the glider may be 
stable or unstable depending upon its initial sidewise 
displacement, and at angler, below this region unstable 
oscillations and divergences are encountered. 

The stability charts indicate a close degree of cor- 
relation between the theoretical and the flight results. 
The dihedral range of the various stability regions indi- 
cated by theory is corroborated by the flight tests and 
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stability reversal in flight occurs at appr fc&cimat ely the 
same effective dihedral as that indicated "by the stabil- 
ity "boundaries. 

is Tyoical variations of the period of the periodic 

l ^ mode and the damping characteristics of the periodic and 

H aperiodic modeo are presented in figures 9, 10, and 11. 

These figures further amplify the conclusions drawn from 
the boundary charts and indicate that a, high degree of 
stability ip obtained at small positive dihedral angles. 
Ac? shown "by figure 9, maximum damping of the oscillatory 
mode exists for these dihedral angles. Any increase or 
decrease of dihedral away from these small -oositive val- 
ues results in lighter damping of the oscillatory mode 
and eventually results in oscillatory instability. 

Figure 1C shows that an increase in the period of 
the lateral oscillations occurs with reduction of effec- 
tive dihedral angle, resulting in long-period oscillations 
for small negative dihedral angles. Th^se long-period 
oscillations occur in an unstable oscillatory region and, 
hence, are likely to have the appearance of truly diver- 
gent modes that exist at more negative dihedral angles. 
These characteristics may account for the divergence rat- 
ings given the flight test made with small negative ef- 
fective dihedral angles. 

The point of peak stability noted for low positive 
dihedral angles in the theoretical analysis is in agree- 
ment with the flight test data in reference 1, in which 
the steadiest and most stable flights occurred at small 
positive: dihedral angles. As predicted by the theoretical 
analysis, divergent conditions were encountered in the 
flight tests for negative values of effective dihedral. 
Reasonably close quantitative agreement was obtained for 
the period of the oscillatory mode as shown in figure 10, 
and both investigations indicated an increase in period 
with decrease in effective dihedral angle. 



Effect of Sidewise Displacement 

Figures 9 and 11 illustrate the adverse effect of 
sidewise displacement upon the stability characteristics 
of the glider. A decrease in the damping of all stability 
modes accompanies sidewise displacement. Although this 
effect is small for small positive dihedral angles, it 
becomes the predominant factor in the determination of 



22 



the stability for large positive dihedrals and eventually 
causes oscillatory instability or divergent flight. These 
results explain certain phenomena previously encountered 
in the flight tests reported in reference 1 for flights 
with large dihedral angles. During the tests the model 
would fly quite satisfactorily for some period of time at 
large dihedral angles "before unstable oscillations devel- 
oped. 

In order to obtain more flight data on the effect of 
sidewise displacement ^ additional flight tests have since 
been made in which a model equipped with remote controls 
wa? attached to a tunnel support by means of a dyadic 
towline system. This model, although larger and heavier 
than the original test model, portrayed the same essential 
variation of stability with dihedral angle as the original 
model. The model was then flown at level tow and aileron' 
control disturbances were introduced during the flight. 
In these flights the model would fly satisfactorily in 
the center position and would even damp out light control 
disturbances. For large control disturbances involving 
large sidewise displacements, however, unstable oscilla- 
tions would occur. These tests confirmed the theoretical 
premise that the stability of the glider at high values 
of effective dihedral angle is largely dependent upon ini- 
tial sidewise displacement. The results indicate that 
gusty atmospheric conditions may induce instability in 
flights when the dihedral angle of the glider is large; 
whereas complete stability would exist for smoother air- 
flow conditions. The advantages of small dihedrals are 
consequently again emphasized. 

Effect of Towline Length 

The effect of towline length upon the period and 
damping characteristics of the glider is illustrated in 
figures 12, 13, and 14. These figures indicate that in- 
creasing the length of the tovlines seriously lessened 
the degree of damping of the oscillatory and aperiodic 
modes. A reduction of stability was evidenced for all 
dihedral angles for lengthened towlines, and the peak 
damping of the oscillatory mode for 3-span towlines was 
approximately 50 percent of that obtained with 1-span tow- 
lines. The damping of the aperiodic mode was also seri- 
ously lessened by increase in length of the towlines, 
particularly for the higher value of lift coefficient. 
This effect -probably accounts for the discrepancy between 
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the flight results and the calculated result* for a lift 
coefficient of 0.75 and S-spaa towlines. (See fig. 7.) 
Apparently, the damping of the aperiodic mode for these 
flights vas so low, even though positive, that sustained 

> flights could not "be made. 

to 

w In general, increasing the towline length resulted in 

a slight increase in the period of the lateral oscilla- 
tions , 

The flight tests of reference 1 were in good agree- 
ment witn the theoretical study with regard to the effect 
of towline length. The lessening of the glider stability 
and the more or less constant nature of the period of the 
lateral oscillations with increased towline length were 
noted throughout the flight tests. 



Effect of Lift Coefficient 

The theoretical investigation showed that increasing 
the lift coefficient from 0.30 to C.75 slightly increased 
the period of the oscillatory mode and reduced the damp- 
ing of both the oscillatory and aperiodic modes. This ef- 
fect of increased lift coefficient upon the period of the 
lateral oscillations is illustrated in figure 14 and upon 
the damping characteristics is shown in figures 15 and 16. 

Although for 1-span towlines only moderate effects 
of lift coefficients were evident, for longer towlines 
the adverse effects of high values of lift coefficient he- 
came more serious. Because of the lo^er stability asso- 
ciated with the larger value of lift coefficient, it «ra«i 
indicated that flights at large positive dihedrals would 
"become unstable with side^iee displacements smaller than 
those displacements inducing instability at the lower val- 
ue of lift coefficient. This trend was accentuated with 
increased towline length until, with 3~span towline?, 
complete instability occurred at a value of 14 affective 
dihedral for all values of sidewise displacement, as shown 
i n figure 7 , 

The flight data were in good agreement with the theo- 
retical results. The reduction in stability with increase 
of lift coefficient in the flight tests was slight for 
1-span towlines but became noticeable and serious for 
longer towlines. 
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Effect of Lift-Drag Ratio 

The result? of the theoretical survey as shown in 
figure 17 revealed a decided advantage of low lift-drag 
ratios with respect to the damping of the oscillatory 
node. These data apply only to zero sldewi se displace- 
ment "but similar trends occurred for other values of dis- 
placement . 

The results indicate that greatest stability will he 
obtained at low values of lift-drag ratio. The damping 
of lateral oscillations reached a maximum at a lift-drag 
ratio of 4 for the glider investigated. Increasing this 
ratio to a higher va.lu^ resulted in less stability for all 
dihedral angles although as would "be expected, this effect 
was least pronounced for small positive values of dihedral 
angl e . 

The highly advantageous effect of spoilers, which re- 
duce lift-drag ratio, as reported in the flight tests of 
reference 1 may he seen to confirm these data. 



CONCLUSIONS 



A theoretical study of the lateral stability of a 
glider towed by twin parallel towlines indicates the fol- 
lowing conclusions i 

1) A dyadic system of parallel towlines that imposes 
a restraint in yawing will provide satisfactory inherent 
stability for a pilotless towed glider. 

2. The lateral stability will be chiefly influenced 
by the dihedral angle of the wing. A maximum of stability 
will exist over a range of small positive values of ef- 
fective dihedral. Any variation of dihedral from this 
range will lead to progressively less stability and will 
eventually result in unstable oscillations or divergences. 

3. Ghisty atmospheric conditions, which result in 
large sidewise displacements of the glider relative to 
the tug, will reduce the stability inherent in the tow 
system. This effect will be greatest for large positive 
values of effective dihedral. 



4. Long tovlir.es will decrease the lateral stability 
of the glider although stable satisfactory flight will 
still be obtainable for towlines up to three spans. 

5. Increasing the lift coefficient of the glider 
will reduce the lateral stability for all conditions. 
This effect will be slight for short t owlines but will 
increase in magnitude with long t owlines. 

6. Low lift-drag ratios will increase the lateral 
stability of the glider. 



Langley Memorial Aeronautical Laboratory, 

National Advisory Committee for Aeronaut toe 9 
Langley Field, Va* 
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8 Figure 8.- Damping characteristics of the aperiodic mode for a glider towed "by a dyadic tov/line 13 
system. 2 -span towlines. 
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•Typical variation of period of oscillatory mode vith dihedral and sidewise 

displacement. Cl = 0.30, 1 r-span fcowlinea. 
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Typical variation of lamp in- of aperiodic ^odoAg -ith dihedral and sidev/ise displacement. 1 
KL-.30. 2- 8p°oi towlinss. h 



L-372 




^1 

Figure 12.- Effect of towlength on the damping of the periodic modes ^3, of a glider towed by a o*j 
dyadic tow system. C L = 0.?5« 
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Figure 14.- Effect of tovvlongth on the period of the oscillatory mode of a glider towed "by- 
dyadic tow system, y' = 0. 
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